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Abstract 

It is shown that the current-perpendicular-to-plane giant magneto-resistance 
(CPP-GMR) oscillations, in the ballistic regime, are strongly correlated with 
those of the exchange coupling (J). Both the GMR and J are treated on 
equal footing within a rigorously solvable tight-binding single-band model. 
The strong correlation consists in sharing asymptotically the same period, 
determined by the spacer Fermi surface, and oscillating with varying spacer 
thickness predominantly in opposite phases. 
PACS numbers: 72.15.Gd, 73.14.Jn,75.70.Fr 

The oscillatory behaviour of many physical phenomena of magnetic multilayer systems 
manifests itself in the most spectacular way as a function of spacer thickness, but the mag- 
netic layer thickness is relevant^^, too. The most widely studied oscillatory phenomena 
are those connected with either the exchange coupling (J) or the so called giant magneto- 
resistance (GMR) (see Refs. 3 and 4 for a review of the currect understanding of these phe- 
nomena). The exchange coupling is of quantum nature and is well understood in terms of 
such theoretical approaches like: RKKY-type theory! , quantum well states!, tight-binding 
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dei, and free-electron-like oneB. From these approaches as well as experimental result^l 
and ab initio band structure calculations!'!, consensus emerges on that the oscillation peri- 
ods of J are determined by certain extremal spanning vectors of the spacer Fermi surface. 
As regards the GMR, according to the two-spins channels model, one expects a strong influ- 
ence of the exchange coupling (responsible for the mutual orientation of the magnetization 
of ferromagnetic slabs) on the resistivity. The anticipated trend would be to relate the an- 
tiparallel (parallel) orientation with maxima (minima) of GMR. The GMR can be easily 
measured if the relative spontaneous orientation of the magnetizations of the magnetic slabs 
is antiparallel (negative J), since then simply GMR= (-R(O) — R{H))/R{0), where H is the 
magnetic field necessary to switch to the parallel orientation; but GMR remains well defined 
in the opposite case, too. While the latter case makes no problem for a theoretical treat- 
ment, it requires pretty sophisticated handling (atomic engineering) in order to stabilize the 
antiparallel orientation by pinning one of the ferromagnetic slab magnetizationdl^. 

Although the GMR in general is not of quantum origin and contains some ingredients 
which are hard to control (defects, impurities, surface and interface roughness etc.), there 
is one contribution, due to reflections of electrons from quantum well barriers, which is 
of the same origin as the exchange coupling. This quantum contribution has been stud- 
ied and shown to be quite substantial both by first principles computationsS and model 
calculations0'Ei. The aim of the present paper is to confront the GMR oscillations in the 
ballistic regimeSHEl, where only the quantum contribution appears, with those of the ex- 
change coupling. Both quantities are treated on equal footing without any approximations, 
by precise numerical computations. 

It is interesting that in spite of plenty of publications on GMR and J, there has been, 
to our knowledge, only one theoretical attempt^ devoted to a detailed comparison of both 
quantities. The authors of Ref. 7 have used a free-electron-like model, compared the J 
behaviour with that of the current-in-plane (GIF) GMR, and found that the CIF-GMR 
assumes maxima for the parallel orientation. That study could not treat, however, the 
relevant quantities on equal footing since some inconsistency was unavoidable as a result of 



taking into account electron scattering by impurities. As we will see below, in our model 
the maxima, for the current-perpendicular-to-plane (CPP) GMR arise for the antiparallel 
orientation 

In the present letter we put emphasis on selecting the above-mentioned CPP-GMR com- 
ponent of purely quantum origin, which can be directly compared with the exchange cou- 
pling, and in principle, measured in the ballistic regime^. We adopt the rigorously solvable 
tight-binding single-band model Hamiltonian 

^ = E%4^/.^ + EK(?)cU.,, (1) 

with t being the nearest-neighbour hopping integral (|t| is the energy unit) and Vaii) - the 
spin-dependent atomic potential. The systems under consideration are trilayers of the type 
UfF/nsS/rifF , where Uf {ug) stands for the number of the ferromagnetic (non-magnetic 
spacer) monolayers in the perpendicular z -direction. The ferromagnetic slabs are mag- 
netized either parallel or antiparallel to each other. Since the systems are infinite in the 
(x, ?/)-plane and the potentials are only z-dependent, the eigenvalues of the Hamiltonian 
(15) for simple cubic lattices (with the lattice constant = 1) are simply: 

— * 

£'(/c||, r) = e_L(r) — 2(cos + cos ky) , (2) 

where e± are the eigen- values, labeled by r, of a tridigonal matrix of rank 2nf + ns (with 
free boundary conditions in the z-direction). The exchange splittings of the ferromagnetic 
films have been introduced by taking spin- dependent atomic potentials V^- outside the spacer. 
Additionally we have assumed a perfect matching of the minority bands in the whole system 
by putting V| = 0. The exchange coupling is then calculated in the following wayS'S 

= Y.[E{h^) - Ef] ■ e{EF - E{k\\,T)), (3) 

fe||,r 

j = f]n + ^i_^T_^T^ (4) 
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i.e. J is the difference of the grand-canonical thermodynamic potentials for antiparallely 



accurately by means of the special k-points methodEH. 

The way we compute the CPP-GMR is the same exact one as in our recent paper§. 
It uses the Kubo formula and applies an accurate recursion Green's function method to 
trilayer systems sandwiched between semi-infinite ideal lead wires. The method, based on 
Refs. 18,19 and modified by us in Ref. 13, is rigorous. The GMR is defined by 



where T is the conductance, and the superscripts and subscripts indicate the relative 
orientation of the magnetization, and the electron spin, respectively. 

In Fig. 1 we present some typical plots of GMR vs. Ug for different magnetic slab thick- 
nesses (rif). What is easy to see is that for small Uf the curves have a regular quasiperiodic 
behaviour which seems to be disturbed for nf greater than, say, 6. Before we explain this 
puzzling situation let us briefly remind that in the asymptotic limit of large Ug one expects 
a quasiperiodic behaviour of the exchange coupling with a period depending exclusively on 
the Fermi energy which (for a given lattice structure) determines extremal spanning vec- 
tors of the spacer Fermi surface. Fig. 2 examplifies this for the 3-F/105'/3-F system with 
Ep = 2.5. The Fermi surface has the form of e_L-constant contours with occupied electronic 
states inside. The arrow indicates the spanning vector Q with the following components 
[k^ = 71, ky = Q = 0.774, e_L = 1.929]. We quote the numerical values to show that already 
for a relatively small system with 16 monolayers in the 2;-direction e± is close to the asymp- 
totic value 2, i.e. kz = vr. Thus, asymptotically, for Ug — > oo, when the equivalence of all the 
directions is restored and k^ becomes a good quantum number, one gets spanning vectors 
[77,(5, vr] ^-iid, by symmetry, [11,71, Q], i.e. the oscillation period A = tt/Q ~ 4. Incidently, it 
is very easy to predict a period length in this limit, by minimizing with respect to k^ and 
ky the following function: 




GMR 




-1, 



(5) 




(6) 
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which gives Q = A;/^*'- for {h,ky) = (0,0), (±7r,0), (0,±7r), (±7r,±7r). Another fact 
worth mentioning is that the influence of the magnetic slab thickness on the J oscillations 
only leads to some phase shifts without actually changing the period lengthiS. 

In Fig. 3 we show both GMR and J for n/ = 3 and V| = —1.8 and find pretty well 
correlated oscillations with the same period consistent with the strictly calculated Fermi 
surface and the sketchy estimations above. It turns out that this conclusion holds also for 
the other curves in Fig. 1, but to see it one has to take a closer eye at them and allow for 
greater values of to select the asymptotic trend. That has been made clear in Figs. 4 
and 5 for n/ = 7 and 5, with a period of about 10 in the latter case. It is easily seen 
that the GMR does share with the exchange coupling the long period of oscillations but 
has predominantly an opposite phas^, in the sense that for negative (positive) J it takes 
larger (smaller) values than its aymptotic value. This coincidence is hardly visible for small 
Us until the asymptotic behaviour develops, and is partially obscured by the superposition 
of some short period oscillations of GMR of non- RKKY natureEili. 

In conclusion, we have carried out numerical studies of a rigorously solvable tight-binding 
single-band model - treating both the CPP-GMR and the exchange coupling J on equal 
footing - and found that asymptotically both quantities share the same oscillation period 
(consistent with the Fermi surface topology) and have predominantly opposite phases. Since 
according to Ref. 19 the CPP-GMR is rather non-sensitive to impurities, our statements 
should also apply if additionally impurities are taken into account. 

This work has been carried out under the bilateral project DFG/PAN 436 POL (U.K., 
S.K.) and the KBN grants No. PB.872/P03/96/11 and 2P 03B 165 10 (S.K.). We also 
thank the Poznah, and Regensburg Computer Centres for computing time. 
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Figure Captions 



Fig. 1: CPP-GMR of the systems UfF/usS/nfF (where Uf and Ug are the number of 
ferromagnetic and non- magnetic spacer monolayers, respectively). Majority spin electrons 
have the potentials = —1.8 in the ferromagnets (all other potentials are 0), Ep = 2.5. 

Fig. 2: Fermi "surfaces" of the systems under consideration consist of ej_-constant energy 
contours with occupied states inside, where e± fulfils: e± — 2(cos/i;-r + cos ky) — Ep. The 
arrow indicates the extremal spanning vector Q for the 3F/10S/3F system (in the parallel 
configuration, with Ep = 2.5 and V| = —1.8). 

Fig. 3: GMR (solid line) and exchange coupling (J) (dashed line) vs. the spacer thickness, 
for Uf = 3, Ef = 2.5 and = -1.8. 

Fig. 4: The same as Fig. 3, but with Uf — 7 and for higher values of Ug to reveal the 
asymptotic trend of the oscillations. 

Fig. 5: The same as Fig. 4, but with n/ = 5, Ep = 2.1, and = —2. 
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